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Abstract 

As a contribution to the ongoing controversy about the role of higher-order Kerr effect (HOKE) in laser 
filamentation, we first provide thorough details about the protocol that has been employed to infer the 
HOKE indices from the experiment. Next, we discuss potential sources of artifact in the experimental 
measurements of these terms and show that neither the value of the observed birefringence, nor its inversion, 
nor the intensity at which it is observed, appear to be flawed. Furthermore, we argue that, independently on 
our values, the principle of including HOKE is straightforward. Due to the different temporal and spectral 
dynamics, the respective efficiency of defocusing by the plasma and by the HOKE is expected to depend 
substantially on both incident wavelength and pulse duration. The discussion should therefore focus on 
defining the conditions where each filamentation regime dominates. 

PACS numbers: 42.65.Jx Beam trapping, self focusing and defocusing, self-phase modulation; 42.65.Tg Optical 
solitons; 78.20.Ci Optical, 37.10.Vz Mechanical effects of light on atoms, molecules, and ions; 42.50.Md Optical 
transient phenomena: quantum beats, photon echo, free-induction decay, dephasings and revivals, optical nutation, 
and self-induced transparency constants 
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I. INTRODUCTION 



While potentially spectacular applications like rainmaking or the control of lightning QlJ— [31] at- 
tract much attention on the filamentation of ultrashort laser pulses, a controversy has recently been 
raised about the physical mechanism at the root of this phenomenon. Filaments are generally de- 
scribed as a dynamic balance between Kerr self-focusing and defocusing by the plasma generated 
at the non-linear focus Q4j— lZD - The recent measurement of higher-order Kerr effect (HOKE) indices 
of alternate signs in air and argon [8, §| led us to propose that these terms provide the main defo- 
cusing effect, so that ionization and self-guiding are almost decoupled II 1011 . at least in the visible 



and infrared regions ill ill . Based on numerical work, and in the lack of knowledge of their va 
these terms had already been assumed to contribute to defocusin g, bu t only marginally 1112 - 



ues, 



Our unexpected prediction has therefore been actively challenged DJ/7 



16] 



12D. 



The controversy simultaneously focuses on two questions: the validity of the experimental 
measurement of the HOKE indices, and the validity of a filamentation model based on them. It is 
fed by the difficulty to perform quantitative measurements in filaments, due to the high intensity 
within them. This difficulty prevents one to directly test the contribution of the HOKE to fila- 
mentation. In this paper, we address these two aspects, with the aim of making the controversy 
as factual as possible by summarizing the facts and opening questions on this subject. In a first 
section, we establish the methodology used for extracting the HOKE terms from the experiment 
and discuss several potential artifacts in the experimental measurement of the HOKE terms, the 
values of which are critical to evaluate their contribution to filamentation. In a second section, 
we discuss the relevance and the physical consequences of the introduction of the HOKE in the 
description of filamentation. We suggest that the contribution of the HOKE on the filamentation 
process strongly depends on the incident wavelength and pulse duration. More specifically, longer 
wavelengths and shorter pulses are more sensitive to the HOKE, while defocusing by the plasma 
is favored in the case of shorter wavelengths and long pulses. 



II. ON THE MEASUREMENT OF THE HOKE INDICES 



IS 



The key result reported by Loriot et al. [8|,|9D is the saturation and inversion of the instantaneous 
(i.e., at least, much shorter than the experimental resolution of ~100 fs) non-linear refractive index 
at high intensities, which we phenomenologically described as HOKE terms from n^I 2 to ng/ 4 in 
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air (resp., to n w I 5 in argon). These terms have been obtained by a numerical fit on the experimen- 
tal data. The implications to filamentation rely on two aspects of the measurement. On one hand, 
the fact that the Kerr effect can saturate and even become negative at high intensities is necessary 
to enable all-Kerr driven self-guiding as described in I lojl . On the other hand, this result can only 
have practical implications on filamentation if this inversion occurs at an intensity lA nKen =o below 
or close to the clamping intensity predicted by the usual filamentation model relying on the bal- 
ance between Kerr self-focusing and defocusing by the free electrons, i.e., J cr i t ~5x 10 13 W/cm 2 



in air [2^,|21|]. If lAn Ksn =o is higher, the Kerr inversion will not occur in laser filaments, and will 
therefore play a negligible role in filamentation. We shall therefore describe in this section the 
experimental protocol that has been used in order to determine the intensity and extract the HOKE 
indices from the experimental sets of data, but also discuss potential sources of artifact affecting 
the measurement of these negative HOKE indices. 

A. Measurement protocol 

In the experiment reported in [8, 2I, the transient birefringence of a molecular or atomic gas 
sample has been measured using a time-resolved polarization technique depicted in Fig.Q] The 
setup allows to carry out two types of detection: homodyne and heterodyne. In the former, the 
signal is related to the squared amplitude of the birefringence (amplitude- sensitive detection), 
whereas in the latter it is related to the birefringence itself (amplitude and phase sensitive de- 
tection). The heterodyne detection therefore provides the sign of the birefringence. In practice 
it is implemented by introducing between the cell and the analyzer a phase plate producing a 
static birefringence (see [SO). The measurements are based on the comparison between two op- 
tical Kerr contributions induced by the field; the electronic Kerr contribution resulting from the 
deformation of the electronic cloud and the reorientation of the molecular dipole due to molecu- 



lar alignment ||22|1 . respectively. This section provides details about the procedure that has been 



followed in order to extract the HOKE indices from this experiment. 

1. Intensity calibration 

A special attention has been paid to the estimation of the laser intensity experienced by the 
molecules or atoms present within the interaction volume. This intensity was inferred from the 
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Figure 1: Experimental setup for time-resolved birefringence measurements. M: Mirror, BS: Beam Splitter, 
L: Lens (f=20cm), CC: Corner Cube, P: Polarizer, A: Analyzer, S: Beam stop, PM: Photomultiplier. The 
relative polarizations of the pump (PI), probe (P2), and signal-field (A) are shown in the inset. 



measurement of the field-free alignment Ii23l l24h. The last, also named post-pulse alignment, is de- 
scribed by (cos 2 9) (t) — 1/3, with 9 the angle between the molecule axis and the field direction and 
where ( } denotes the expectation value averaged over the thermal distribution of molecules [25]. 
It is well established that both the structural shape of the alignment revivals and the permanent 
alignment are very sensitive to the initially applied laser intensity Ii26ll . Below saturation of the 
alignmentj_the value of (cos 2 9) at the revivals of alignment increases linearly with the applied 
intensity 11251 12711. whereas between revivals (i.e., for permanent alignment) it grows first quadrat- 
ically and then linearly with the intensity. Field-free alignment therefore provides an accurate and 
unambiguous estimation of the laser intensity in the gas. 

Figure[2] shows the time-resolved birefringence signal of 2 at different laser energies recorded 
with an heterodyne detection. Pure heterodyne detection provides a post-pulse signal propor- 
tional to the convolution of the probe intensity I pr (t) with ((cos 2 9}(t) — 1/3) [8|, 



28] 



SU(*) oc J pr (t) <8> ((cos 2 9}(t) - I) . 

4n n e n V 6/ 



(1) 



4ra e V 3 

with Aa the polarizability anisotropy, p the gas density, no the linear refractive index, e the 
dielectric constant of vacuum, and ® denotes the convolution. The permanent alignment offsets 
the baseline for positive delays. This offset increases with the intensity. Because of the intensity 
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Figure 2: Pump-brobe birefringence signals (black solid lines) recorded in static cell filled with 0.1 bar of 
O2 at room temperature. The energy of the pump pulse is 26 (a), 56 (b) , and 96 [iJ (c). Numerical fits of the 
postpulse signals (red solid lines) used to evaluate the effective intensity in the experiment: 7=8 (a), 20 (b), 
and 34 TW/cm 2 (c), respectively (see text). The simulations are based on Eq.Q] The horizontal dashed-lines 
represent the zero signal. 

profile of the pump and probe beams, the signal measured in the experiment results from a spatial 
averaging. Using a space-averaged calculation (i.e., a 3D model) we have checked numerically 
that the volume effect can be adequately taken into account in a simpler ID calculation just by 
using an effective intensity. In fact, considering a gaussian beam profile, a peak laser intensity 
/peak, and a crossing angle of about 4° between the pump and probe beams [8, 2I, the field-free 
alignment signal integrated over the volume can be approximated by the signal produced at the 
effective intensity I defined as 



'peak 



/1.7. 



(2) 



This approximation allows to save computer time when fitting the experimental data. 

Figure[2]also shows the results of the simulations that have been used to fit the effective inten- 
sities. The temporal envelope I pi (t) of the probe beam has been described by a gaussian function 
of duration slightly above 100 fs (FWHM) in order to account for the crossing of the two beams. 
The fact that the simulations reproduce very well both the revivals and the permanent alignment 
supports the analysis based on the effective intensity. 
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Figure 3: Effective intensity / fitted on the postpulse alignment signal in O2 (red solid squares) compared to 
the intensity estimated from the measurement of the pulse/beam characteristics (black solid circles). Linear 
regression of the fitted intensity on the measured pulse energy (red solid line). 

The intensity determined as described above has been double-checked by comparing its value 
to the one estimated from the measurement of the beam waist, the incident energy, and the pulse 
duration. The estimation was supported by i) the limited gas pressure in the static cell (0.1 bar, 
except for measurements at or below 1 TW/cm 2 ), ii) the limited power of the incident beam 
P = 1.8GW<C P cr = 80 GW, where P cr is the critical power for nitrogen, as well as the low 
value of the nonlinear phase accumulated during propagation (less than 0.1 and 0.14 rad for N 2 
and 2 , respectively) and (iii) the focused geometry (J = 20 cm), which limits self-channeling that 
would induce changes in the beam profile. This ensures that the propagation was mostly linear in 
the experiment, allowing to estimate the intensity therein from the measurement of the beam waist 
conducted at low energy. Figure[3]reports the effective intensities obtained by the two independent 
methods in the case of nitrogen. The error bars on the measured intensity (solid circles) results 
from the uncertainties on the measured energy (±5 %), pulse duration (±10 %), and beam waist 
(±15 %). The red solid line corresponds to a linear regression of the intensities fitted on the post- 
pulse alignment signal. It allows to estimate an uncertainty of ±10 % (at 3 standard deviations) in 
the determination of the intensities. We can therefore exclude flaws in the intensity determination 
beyond this order of magnitude. 
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2. Determination of the lowest-order Kerr index ( n 2 ) 



The measurement of the lowest-order Kerr index has been conducted at low energy in order to 
avoid any influence from the HOKE terms. In order to ensure a good signal-to-noise ratio despite 
the weakness of the birefringence signal, an homodyne detection, that offers a better sensitivity 
than the heterodyne one, has been employed. In the case of a homodyne detection, the signal is 
given by BSD 

/ 3pAa / 2 . . 1 



\4n e V 3, 
+ j^2 r ° ss A*)) , (3) 

where the lowest-order electronic Kerr response n 2 0SS has been added to the retarded rotational 
response resulting from the alignment. Here, it is worth mentioning that the signal measured in 
our experiment results from the cross-coupling between two distinguishable laser beams, namely 
the pump and probe beams. The coefficients ?if oss = 2 x n 2 , with n 2 the standard self-induced 
Kerr index, therefore describes the non-linear refractive index due to cross-Kerr effect. The corre- 
spondence between cross-Kerr and Kerr indices is given in \^\. Finally, the factor 2/3 results from 
the different values of the Kerr index experienced by the probe field along directions parallel (||) 
and perpendicular (_L) to the pump field, respectively, with n 2 ,, = 3 x n 2± . This relation is valid 
when the intrapulse alignment can be neglected so that the medium can be viewed as isotropic 



during the pump excitation B29I1 . The approximation is justified for the investigated molecules and 
the relatively short pulse duration used in the experiment compared to the rotational period. For 
instance, in N 2 or O2, the orientational Kerr contribution to n 2 calculated from the elements of the 



hyper-polarizability tensor 13011 is less than 5% at the maximum peak intensity investigated in the 
experiment. 

Since the permanent alignment can be neglected in the low intensity regime, both rotational 
and electronic contributions exhibit a linear dependency on the applied pump intensity 7, which 
therefore acts as an amplitude factor on the overall signal. n c 2 mss was hence determined indepen- 
dently from the knowledge of this intensity. Figured shows the time-resolved birefringence signal 
recorded in O2 at low intensity. For all recorded data, the value of n 2 has been evaluated through 
deconvolution of the rotational response from the signal. First, the postpulse signal was adjusted 
by Eq.[3]with nf oss = 0, using a low intensity value as a fixed parameter and an amplitude factor 
as a free parameter. This allows to evaluate the contribution of the rotational term close to the zero 
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Figure 4: Homodyne birefringence signal (black solid squares) recorded in 1 bar of O2 at room temperature 
and numerical fit (red solid line). The effective intensity of the pump is 1 TW/cm 2 . In the inset, the total 
birefringence (red solid line) resulting from the electronic (blue dotted line) and rotational (green dashed 
line) Ken - contributions. 

delay, and only relies on the knowledge of the molecular polaris ability (A a = 4.6 and 7.25 a.u. 
for N 2 and Oo, respectively). Next, the signal was adjusted using n| ross as a free parameter. As 
described in Q8|, lSQ , measurements in argon have been calibrated with the postpulse signal recorded 
independently in N 2 in the same experimental conditions. The data presented in Table 1 of Q8L 
result from statistics performed over 50, 30, and 8 data samples recorded at low intensity (i.e., 
I < 1 TW/cm 2 ) in N 2 , 2 , and Ar, respectively. 

3. Determination of the higher-order Kerr indices 

All HOKE indices have been determined by the use of an heterodyne detection that is phase 
sensitive and therefore allows to recover the sign of each term. When considering HOKE terms 
up to the fifth power of the applied intensity /, the heterodyne birefringence signal is given by [8J] 

State (*) oc J pr (t) ® (An wt {t) + jV 2 ross /(t) 
4 fi 



9 8 W 11 



n cros S/ ^4 + lP_ ra cro SS/ ^5 j ? (4) 
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n 



2j 
C 2j 



" i 



1.04 



1.30 



1.45 



10 



2.66 



TABLE I: Correction factors for the volume effect defined as c%j 
HOKE terms and its corresponding effective value, respectively. 



n, 



2j 



,ss /n" oss , with n| r ° ss and raf ° ss the 



'2j 



'2j 



2j 



with An rot (t) = 3pAa/4n e ((cos 2 6*)(t) — |) . We have generalized the relation 712*, = 
(2j + 1) x n 2 j x with j G N*, verified for n 2 and n 4 Bill , to higher orders. As mentioned in the 
previous section, this approximation results from neglecting the intrapulse alignment considering 
the medium as isotropic during the interaction with the pump. 

Although the determination of nf oss is straightforward, since it is independent from /, the 
evaluation of the HOKE indices is complicated by the fact that spatial averaging depends on the 
non-linearity order. As in Sec. lII A ll in order to avoid the prohibitively large use of computer time, 
the deconvolution of the HOKE from the signal has been achieved by fitting the experimental data 
with a ID model. However, in order to evaluate the influence of the volume effect, different 3D 
simulations of Eq.@]have been preliminary performed. First, Eq.@]has been truncated to the second 
power of I (i.e., n.g mss , • • ■ ,nfQ SS = 0) and then spatially integrated. Second, the ID model, in 
which ri4 ross was replaced beforehand by the effective index r^ ross , has been used to fit the previous 
numerical result using the effective intensity defined in Sec. lII A fl as a fixed parameter and n°™ ss 
as a free parameter. For the next HOKE index rig ross , the same procedure has been applied. The 
result of the 3D calculation including n c ™ ss and ?7,g ross has been fitted with the ID model with 
^cross g xe( j 5 f^ ross being the free parameter. This approach has been repeated successively for each 
HOKE index up to nf oss . In order to check that the ratios between the HOKE indices and their 
respective effective values were independent from the intensity, different numerical tests have 
been performed over the intensity range considered experimentally. The ratio between the HOKE 
indices and the effective values that account for the volume effect are given in Table|IJ Since both 
n 2 and the alignment depend linearly on the intensity, the correction factor for n 2 is c 2 = 1 (i.e., 

,4 2 — n 2 '■ 

FigureO displays the time-resolved birefringence signals recorded in N 2 from low to high in- 
tensity. The simulations have been performed using the n|™ ss coefficients first reported in Q8|, [9Q , 
determined as detailed above. Like in 2 and Ar, the electronic Kerr contribution to the birefrin- 
gence signal (i) scales linearly with the intensity at low intensity where n| ross dominates, (ii) then 
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Figure 5: Intensity dependence of the birefringence signals (black solid lines) recorded in static cell filled 
with 0.1 bar of N2 at room temperature. Numerical simulations of Eq.@](red solid lines) where nfj° ss and I 
have been replaced by the effective values r£[° ss and I, respectively (see text). 7=13 (a), 20 (b), 27 (c), 35 
(d), 39 (e), 42 (f), 45 (g), and 49 TW/cm 2 (h). 

saturates at moderate intensity, and finally (iii) reverses its sign and becomes highly non-linear at 
high intensity where the HOKE dominates. To avoid correlation between the n|™ ss coefficients in 
the least square fitting procedure orders were fitted successfully. rf™ ss was first adjusted with nf 055 
fixed and n^j°> 4 = 0, from a set of 14 (resp., 10 and 11) data recorded in N 2 (resp., 2 and Ar) at 
an effective intensity I ^27 TW/cm 2 (resp., 20 and 24 TW/cm 2 ). Due to the predominant value of 
7If oss compared to n^ oss , it has not been possible to isolate an intensity window where these two 
indices could be fitted independently. They have therefore been determined simultaneously, with 
n cross an( j ^cross fj Xec j anc j ^cross = o, from a set of 24 (resp., 10 and 26) data recorded in N 2 (resp., 
2 and Ar) at 27 < 7 < 50 TW/cm 2 (resp., 20 < 7 < 35 TW/cm 2 and 24 < 7 < 45 TW/cm 2 ). For 
the same reason, the last term r^ r oss , necessary only in argon, has been fitted together with n^° ss 
andrf c 8 ross . 

In order to reduce the errors bars, the records mentioned above have also been used in order 
to fit the HOKE indices on the two-dimensional data displaying birefringence signal as a function 
of intensity and time delay. Figure[6] displays a comparison between the experimental data set 
and the numerical simulations resulting form the fitting procedure. The good agreement between 
observations and calculations confirm the values of the HOKE indices retrieved as detailed above, 
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as evidenced in particular by the fact that they overcome the rotational response (see positive signal 
at positive delays for N 2 and 2 ) at large intensity. 




Figure 6: Time-resolved birefringence signals recorded for different effective intensities I in 2 (a), N 2 
(c), and Ar (e). Numerical fits of the 2 (b), N2 (d), and Ar (f) data normalized for each intensity to the 
maximum amplitude of the postpulse signal (not shown), except for argon. 



B. Is the observed inversion due to negative HOKE? 

Negative non-linear indices generated by high-power laser pulses are generally attributed to a 
contribution of the free electrons, given by the Drude model Q5[| : 

^Tlplasma ~Z j (5) 

AriQPcrit 

2 

where p is the electronic density, n is the linear refractive index, p crit = eon ^ w i [ s the critical 
plasma density (p cr i t = 1.75xl0 27 m -3 at 800 nm), e is the permittivity of vacuum, m e is the 
electron mass, lui is the laser angular frequency, and e the elementary charge. 

As described in the previous section, our measurements were performed in a pump-probe con- 
figuration where the pump beam intensity is sufficient to partially ionize the gas, so that the free 
electrons contribute negatively to the refractive index. To cancel this free electron contribution, 
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we measured the transient variation of the birefringence rather than the variation of the refrac- 
tive index itself. This approach is supported by the general belief that non relativistic plasmas 
are not birefringent, even at the time scale of the laser-medium interaction Ii29ll . However, since 
this statement has never been proved rigourously to date, we discuss below the potential plasma 
contribution on each time scale. 



1. Postpulse plasma contribution 



The birefringence induced by a high-intensity ultrashort pulse in argon is not maintained after 
the pulse has passed B2M34I1. A direct contribution of the plasma to the medium birefringence 
would have the same lifetime as the plasma, i.e., at least several ps [350. It can therefore be 
excluded. 
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Figure 7: Time-resolved birefringence signals (solid lines) recorded for different effective intensities 7=10 
(a), 21 (b), 28 (c), and 43 TW/cm 2 (d) in argon. Numerical fits are displayed as red dotted lines. 

This negligible birefringence of the plasma is further confirmed by the temporal dependence of 
the birefringence in our experiments. In the case of argon (see Fig. [7]) the birefringence goes down 
to zero for time delays larger than ~200fs, i.e., as soon as the two pulses do not overlap anymore. 
In the case of oxygen and nitrogen (see Figs.[2]and[5l respectively), the birefringence observed 
at positive delays between the revivals of alignment is perfectly reproduced by the molecular 
permanent alignment described in Sect. HI ATI 
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Figure 8: (a) Pulse intensity (red line) and electron density (black dotted line) calculated in argon along 
the pulse duration (100 fs) for a peak intensity of 28 TW/cm 2 . (b) Corresponding refractive index produced 
by a positive Kerr effect (blue dotted line), an hypothetical negative Kerr effect resulting from accumulated 
free electrons (black dashed line), and both contributions (red solid line), (c) Heterodyned signal resulting 
from both contributions. 

2. Intrapulse plasma contribution 

Non-birefringent plasma could however be expected to reduce the birefringence by inducing 
an equal refractive index change on both axes. However, free electrons accumulate over the pulse 
duration, so that this potential contribution to the non-linear birefringence should be asymmetric 
in time, as depicted in Fig.[8l In contrast, for a symmetric pulse like we used in[|81], the Kerr 
contribution to the refractive index is symmetrical in time. The relative contribution of the plasma 
to the birefringence would be therefore defined by the asymmetry of the birefringence profile 
around time t — 0. Such asymmetry is not observed in our experiment, as evidenced for example 
by the data in argon displayed in Fig. [7] 

Besides the free electrons, the plasma is made of ions. Their susceptibility can be estimated 
according to Sprangle et al. [13611 : 



(3) 
A ion 

(3) 
/^neutral 



IP 



neutral 



IPin 



(6) 



where IP a denotes the ionization potential of species x, as detailed in Table [Til Assuming that the 
refractive index variation of the ions is three times larger along the laser polarization axis than 
perpendicular to it, as is the case for the neutrals, the contribution of the ions is about 10% to 
20% lower than that of neutrals for the same partial pressure. At the maximum intensity used our 



experiment, the ionization is restricted to 1% at most j37\\ . so that the ions cannot, in any case, 



contribute to more than 0.2% to the observed birefringence. Moreover, as discussed above for the 
electrons, any contribution from the ions would result in an asymmetric temporal profile. 

Finally, ionization can deplete the ground state population to the benefit of excited bound and 
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Species 


N 2 


o 2 


Ar 


IPneutral(eV) 


15.6 


12.07 


15.8 


IPion(eV) 


27.12 


24.14 


27.6 


x (3) 

(3) 
-^■neutral 


0.19 


0.125 


0.185 



TABLE II: Ionization potentials of neutral and ionized nitrogen, oxygen, and argon. 



continuum states. Both processes can result in a large modification of the refractive index 1138- 
I40L Moreover, for intensities close to the inversion of the refractive index, the ionization takes 
place at the frontier between multiphoton and tunneling ionization regimes. Besides depletion 
of the ground state, ionization might contribute to the HOKE through the fast-moving electrons 
recolliding with the atomic or molecular core, as described by the three-step model 114 ill . The 
resulting oscillation of the induced dipole responsible for high harmonic generation I142ll could lead 
to large nonlinearities in the Kerr effect. Since we consider a process occurring within an optical 
cycle, this potential effect can be considered as instantaneous as compared with the time scale of 
the pulse. Calculations of the dipole induced by a strong laser field could therefore contribute to 
interpret the effect observed in our experiment. 



3. Two-beam coupling 

Another artifact that could be raised is the energy exchange between the two crossing laser 
pulses, known as two-beam coupling [29, 43]]. First, we should mention that the signal provided 



by heterodyne detection used in the high-order Kerr measurement is in principle free from any 
two-beam coupling contribution. In fact, pure heterodyne detection results from the difference 
between two data set recorded in the same conditions, except for opposite phases of a local oscil- 



lator [44, 



4511 . Second, for femtosecond pulses with identical spectra, two-beam coupling requires 
a frequency chirp and a finite time response of the nonlinear refractive index I146ll . The very 
fast excitation time associated with non-resonant excitation combined with the small residual fre- 
quency chirp of our pulses can only lead to a marginal amount of two-beam coupling through the 
purely electronic Kerr response. Finally, and most importantly, the time profile resulting from 
two-beam coupling should be asymmetric with respect to the pump-probe delay. This was not ob- 
served in our experiment, as previously mentioned (See e.g. Fig. [7]). Consequently, the influence 
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of two-beam coupling on our experimental results can be confidently disregarded. 



HI. HOKE- AND PLASMA-DRIVEN FILAMENTATION REGIMES 

Improving models by considering higher-order of the relevant processes is a quite natural ap- 
proach in all branches of physics. Higher-order optical non-linearity is well-known, giving rise 
e.g. to higher harmonics Ii29ll and should therefore be included in the propagation equations, un- 
less the it is proven that it induces a negligible effect in the considered situation. In the context 
of filamentation, such development was already tried several years ago |l2- 16], although the sig- 



nificance of those works was limited by the lack of experimental knowledge of the magnitude of 
these terms. 

This magnitude constitutes the key question arising about their role in laser filamentation. Con- 
sidering the wide range of wavelengths, pulse durations, incident energies and focusings investi- 
gated in the last ten years Q4j— lZD , we expect that no unique dominant effect can be identified for 
all conditions. Rather, four regularizing terms exist in the full non-linear Schrodinger equation 
(NLSE) [Q], with relative contributions depending on the experimental conditions. Three of them 
(defocusing by the HOKE and by the free electrons, as well as diffraction) are spatial, while 
group-velocity dispersion (GVD) is temporal. 

The latter two are independent on intensity. Diffraction cannot be the dominant regularizing 
factor in filamentation, since it can only occur beyond the self-focusing P c , defined as the power 
required for the self-focusing to dominate diffraction. Furthermore, due to the limited bandwidths 
at play and the relatively low dispersion of usual gases in the near infrared, the contribution of 
GVD is very small in usual filamenting conditions, especially on the short distances of the labora- 
tory. However, for few-cycle pulses centered in the UV, the bandwidth and the dispersion increase, 
so that GVD may play a substantial role as it does e.g. in water I47L 

In most usual cases, defocusing by plasma and/or the HOKE will therefore be the main regu- 
larizing terms. The relative contributions of the HOKE and the plasma-induced refractive index 
change can be characterized as the ratio of the corresponding terms in the NLSE, at any location r 



and time t 111 ill : 



^oPcrit 

The overall action of both effects on the whole pulse duration is described by the pulse- 
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Figure 9: (a) Dispersion of the fifth-order Kerr index (714) 14811 and (b) electron density generated in air 
by a 30 fs pulse of constant intensity 50 TW/cm 2 . Note that the vertical scales are respectively linear and 
logarithmic. 



integrated value of £: 



/IE,>2^/(r,^'| x \e(r,t)\dt 



I 



P(r,t) 



2n 0Pc „ W,t)\dt 

where e is the reduced scalar envelope defined such that \e\ 2 = I 111 111 . 



(8) 



Both HOKE and plasma defocusing generate a negative non-linear refractive index. However, 
they differ by three aspects. First, multiphoton ionization requires 8 photons in O2 and 1 1 photons 
in N 2 and Ar, so that the corresponding non-linearity is of 8th- or 11th order. In contrast, over 
the intensity range relevant for filamentation, orders up to 4 (resp. 5) only have to be considered 
in air (resp. argon), as discussed in Sec. lII A3l Secondly, the Kerr effect has a time constant of a 
few fs at most, shorter than the pulse duration, while the plasma density accumulates over whole 
duration of the same pulse. Finally, both terms increase in absolute values for higher frequencies, 
but the Kerr effect shows a relatively slow dispersion over the visible spectral range 1481] . while 
the increase of the ionization rates is much faster [5], as illustrated in Fig.|9]in the case of air. 

One can therefore expect that the plasma will tend to provide the dominant defocusing contri- 
bution on the short- wavelength side, while the HOKE should be favoured on the long-wavelength 
side. Indeed, we observed this transition in a recent numerical work [11]. As a consequence, one 
can consider as a general rule that the regularizing process will be defocusing by the free electrons 
for short wavelengths (A < 300-400 nm) and by the HOKE for longer wavelengths. The discus- 
sion should therefore focus on defining the border between "long" and "short" wavelengths, rather 
than on the existence of the HOKE. Still, the position of this transition depends little on the value 
of the HOKE indices since the ionization coefficients cover several orders of magnitudes over the 
spectrum. 

Similarly, as a consequence of the different temporal dynamics of the plasma and the HOKE, 
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Figure 10: Numerical simulation of the space-time dynamics of the on-axis ratio £(r = 0) of the respec- 
tive contributions of higher-order Kerr terms and plasma to defocusing (see text for detail). The white 
regions correspond to irrelevant values due to the low intensity (I < / max /e 2 ). Pulse durations are (a) 30 fs 
d^ D (r = 0) = 3.3 at z = 98.7 cm) and (b) 250 fs (£ min (r = 0) = 0.5 at z = 99.8 cm). Note the asymmetric 
temporal shape due to the plasma accumulation, contrasting with the symmetric pulse shape related with 
the domination of the instantaneous HOKE in panel (a). 



the pulse duration will strongly impact their respective contribution to defocusing. The longer the 
pulse, the more efficient the accumulation of plasma, and consequently its relative contribution 
to defocusing. Indeed, as shown in Fig.fTOh. numerical calculations using the model described 
in 111 ill and relying on the generalized non-linear Schrodinger equation and ionization rates given 
by the multi-species generalized Keldysh-Perelomov, Popov, Terent'ev (PPT) formulation pfl, the 
plasma contribution is marginal (£ 3> 1) over most of the duration of a short pulse. It only 
plays a significant role at the very tail of the pulse, where the intensity has decreased close to 
zero. Its contribution to the overall pulse propagation is therefore negligible, as evidenced by a 
high on-axis value of i min {r = 0) = 3.3 for a 30 fs pulse. In contrast, for a longer pulse (250 
fs, Fig.flOb). the plasma efficiently accumulates earlier in the pulse and is significant already in 
its high-intensity region. It therefore contributes significantly (£ m j n (r = 0) = 0.5 for a 250 fs 
pulse) to the propagation of long pulses, resulting in particular in an asymmetric temporal pulse 
shape. The published values of the HOKE indices suggest that the border between the two regimes 
lies at a few hundreds of fs at 800 nm. The two regimes are therefore experimentally accessible 
with the current laser technology, which may explain the contradicting results obtained among 
the available experimental work. Furthermore, the transition between them is smooth, so that 
in adequate conditions both processes contribute with similar magnitudes. Again, the discussion 
should focus on the domains of the long- and short pulse regime rather than on seeking for one 
single universally dominating mechanism. 
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Figure 11: Schematic "Phase diagram" of the dominating mechanisms at play in filamentation 

The different regimes discussed above can be displayed schematically as a "phase diagram" 
(Fig.QT]). In this graph, the axis have been left blank on purpose to ensure the generality, inde- 
pendently of the measurement of the HOKE terms. From a practical point of view, the transition 
between the regimes corresponds to the equality of the contributions of plasma and HOKE to defo- 
cusing (£ = 1). The historical trend over the last 10 years includes a substantial shortening of the 
pulses, from 100-200 fs in the late 1990's and early 2000's to ~30 fs nowadays. According to our 
results, this shortening corresponds to a transition between "long" and 'short" pulses, i.e., from 
plasma- to HOKE-regularized filaments, which might explain why HOKE had been considered 
marginal up to now. 

IV. CONCLUSION 

As a conclusion, we have discussed the details of the recent measurement of the HOKE indices 
[8|, |9J] and excluded a range of potential artifacts, confirming the reliability of these experimental 
data. Furthermore, we have shown that, due to different temporal and spectral dynamics of the 
plasma- and HOKE- induced defocusing, the former should be dominant for long pulses at short- 
wavelengths, while the latter will dominate for short pulses at long wavelengths. We therefore 
suggest that the controversy about the mechanism of filamentation should turn from a qualitative 
discussion of which effect dominates, to a more quantitative discussion about the domains of 
relevance of each process. 
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